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A new model for three dimensional distribution of atomic hydrogen gas in the Milky Way is 
derived using the 21cm LAB survey data. The global features of the gas distribution such as spiral 
arms are reproduced. The Galactic plane warps outside the solar orbit and the thickness of the gas 
disk flares outward the Galaxy. It is found that the mass of atomic hydrogen gas within a radius of 
20 kpc is 4.3 x 1O 9 M . 



I. INTRODUCTION 

The distribution of gas in the Milky Way reveals the 
global structure and dynamics of the interstellar medium 
[TJ [2]. Furthermore, the fragmentation of cosmic rays 
during their propagation within Galaxy is caused by in- 
teractions with interstellar gas. Besides, the diffuse 7- 
ray components produced by decays of neutral pions and 
bremsstrahlung of electrons and positrons are correlated 
to the gas distribution in the Galaxy. One of the ma- 
jor constituents of the interstellar gas is atomic hydrogen 
which is traced by its 21cm line emission. Several models 
for distribution of atomic hydrogen in our Galaxy have 
been constructed [3-5 . However, recent all sky 21cm sur- 
veys with high angular resolution, strong sensitivity and 
large velocity range motivate us to devise a new model 
with detailed features. 

In this note, we construct a model for three dimen- 
sional distribution of atomic hydrogen gas in our Galaxy. 
To that end, the Leiden-Argentine-Bonn (LAB) survey 
data [6] is used. This survey merges the Instituto Ar- 
gentino de Radioastronomia (IAR) southern sky survey 
[3(8] with the Leiden/Dwingeloo Survey (LDS) [9]. The 
data are corrected for the stray radiation. This survey 
has an angular resolution of 0.6° and a velocity sampling 
of 1 km/s, with velocity range of (-450,400) km/s. It is 
presently the most sensitive 21cm line survey with the 
most extensive spatial and kinematic coverage. 

The derived distribution of atomic hydrogen gas can be 
applied for estimating the interaction rate of cosmic rays 
during their propagation as well as evaluating the diffuse 
gamma rays emission. Small scale features of the gas 
distribution manifest themselves in diffuse gamma ray 
sky maps. These structures can be traced in high angular 
resolution maps of the Fermi gamma ray telescope. 

This paper is organized as follows; In section [TT| the 
properties of 21cm line and its connection to the atomic 



hydrogen number density is reviewed. In section III the 
method of deriving the gas density in the outer and inner 
part of the solar orbit and at tangent points is explained. 
We also discuss the assumptions on the rotation curves. 



In section [IV] we present the general properties of the gas 
distribution and finally, we conclude in section |V} 



II. THE 21 cm LINE EMISSION 

Hydrogen atom emits line radiation with wavelength 
A = 21.1 cm (v = 1420.4058 MHz) through a hyper- 
fine transition when the spins of electron and proton flip 
from being parallel to antiparallel. Since the transition 
probability is too small, collisions have enough time to es- 
tablish an equilibrium distribution of hydrogen atoms in 
the upper and lower states labelled 2 and 1 respectively. 
Thus applying the Boltzmann distribution the ratio of 
the number of atoms in these states is given by 
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where #2 and gi are, respectively, the statistical weights 
of the upper and lower levels with ratio of #2/^1 = 3. The 
excitation temperature T s is called the spin temperature 
and under most circumstances T s ^> hv^jk = 7x 10 -2 K. 
The radiation transfer equation in terms of the radiation 
intensity I V1 namely the radiant energy per second per 
unit area per steradian per bandwidth, can be expressed 
by (for more details see [T0UT3] ) 
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The increase in intensity, in traversing dx, is k u /4ir where 
k v is the emissivity of the plasma. The decrease in in- 
tensity in the same distance increment is Xvlv where \v 
is the absorption per unit path length. 

The emissivity and absorption per unit frequency in- 
terval in terms of the line width of the neutral hydrogen 
profile 5v are 

k v = n2(x)A 2 ihvo/5v, (3) 
Xv = ^{n 1 B 1 2-n 2 B2i)hvQ/5v. (4) 

47T 

The Einstein's coefficients A 2 \, B12 and .821 are the in- 
trinsic properties of atoms and satisfy the following rela- 
tions 



2 



Using [T] and [5| the absorption Q can be rewritten as 
follows 
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where we have used the fact that the number density 
of hydrogen atoms in the lower state is 1/4 of the total 
hydrogen number density. The optical depth in terms of 
absorption \v is 



III. DERIVATION OF ATOMIC HYDROGEN 
NUMBER DENSITY 

The atomic hydrogen number density at a given helio- 
centric distance r, longitude / and latitude b is derived 



from (10) as 
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where 5r is the path length. The optical depth, itself, is 
used to define the brightness temperature 
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The line width Sv is caused by the motions of hydrogen 
gas with respect to the observer. Thus it can be replaced 
by ^-Sv r where v r is the gas radial velocity relative to the 
Sun. Assuming that the motion of gas around the Galac- 
tic center is purely circular, the gas located at Galacto- 
centric radius R and altitude \z\ has a radial velocity with 
respect to the Sun which is given by 

Vr= fOJSA_^L\ RQS]nlcosb (9) 

V R Rq J 

In the above, 0(R, z) is the rotation velocity of the gas, 
I and b are, respectively, the longitude and latitude. Eq. 
([6| is then written as 



where I = 0, 6 = corresponds to the Galactic center. 

The rotation velocity of the gas 0(R, z) away from the 
Galactic plane becomes smaller than that of the gas at 
the underlying disk. At small Galactocentric radii the al- 
titude dependence of rotation curves is prominent while 
it becomes less important at large values of R [17] . More- 
over, the vertical extension of hydrogen gas is small in- 
side the solar circle and it increases outward the Galaxy. 
For these reasons, it is a valid assumption to ignore the 
lagging rotation. 

Inside the solar circle (R < Rq) we use the rotation 
curve of [18] which is fitted by a polynomial of the form 



6(R) 



A R n 



(14) 



The coefficients A n are obtained by assuming Rq = 
8.5 kpc and = 220 km/s. For R > Rq there is a 
general consensus that it is a fair approximation to as- 
sume a flat rotation curve with = 0q [19] [20]. It is 
worth noting that the angular velocity 
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where the constant C is 
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3hc 3 A 21 (11) 
1.823 x 10 18 cm- 2 R- 1 (km/s)- 1 . 



is always positive and increases toward the Galactic cen- 
ter. Radial velocities giving negative uo are forbidden. 
They correspond to the peculiar motions of the local gas. 
The derivative of radial velocity with respect to helio- 



The column density is obtained by integrating eq.dTob centric distance ^ is computed by using the chain rule 



over path lengths and radial velocities along a line of 
sight. 
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The measurement of the brightness temperature X5 over 
a large range of radial velocities and directions in the sky 
is provided by the combined LAB survey [6]. 

The spin temperature is much greater than the bright- 
ness temperature for optically thin (r v <C 1) 21cm 
line emission. However, toward extragalactic sources T s 
varies from 40 to 300 K, depending on the location and 
velocity [14]. Moreover, outside the solar circle up to 
Galactocentric radius of 25 kpc the spin temperature is 
in the range of 250 to 400 K [15] . We assume a glob- 
ally constant spin temperature equal to 150 K, which is 
also the maximum observed X5 in the LAB survey data. 
Although using different values of T s may change the col- 
umn density of atomic hydrogen [16] , the whole structure 
of the gas remains unchanged. 
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Therefore, the gas number density at every given (r, /, b) 
is obtained by inserting eq.(16) in eq.([l3|. The heliocen- 
tric distance r associated to[I, v r ) or (Z,6, i?) is deter- 
mined by 
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In the inner part of the solar orbit there are two kine- 
matically allowed distances (except for the tangent points 
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R = Rq I sin Z | where they coincide) and for the outer part 
there is only one. Although there is no distance ambigu- 
ity at tangent points, because is zero eq.(13) fails to 
determine the gas density. We describe the method of ob- 
taining the gas density at tangent points in section [Til A 



In the case of distance degeneracy the observed inten- 
sity must be distr ibuted among the near-far points as ex- 
plained in section III B Note that is also zero toward 
the Galactic center / anti-center and right above/below 
the Sun. The number density of points in these directions 
is calculated by linear interpolation between n#(r, b) of 
nearby points. 



A. Tangent Points 

The closest point to the Galactic center at every given 
direction rcosb = Rq cos/, has extreme positive (in the 
first quadrant) or negative (in the fourth quadrant) radial 
velocity which is called terminal velocity v t . Due to veloc- 
ity dispersion cr v , the velocity profiles do not have a sharp 
cutoff at tangent points. The emission from the tangent 
points is ideally a bivariate Gaussian in altitude and ve- 
locity. But emission from the nearby radii are not well 
separated in velocity because of the velocity dispersion. 
Atomic hydrogen gas in the vicinity of the tangent point 
has radial velocity with \v r \ > \v t \ — cr v . The number 
density around the tangent point is obtained by dividing 
the emission from this velocity range (\v r \ > \v t \—a v ) by 
the corresponding path length [3]. In the first quadrant, 
it is 



n H (r,l,b) 



cTsi: 



In 1 



T b (l,b,v r ) 



dv 



r2(v t - °v) - n{v t - a v ) 



(18a) 



and in the fourth quadrant it is 
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In the above, r\ and r2 are, respectively, the near and 
far heliocentric distances associated to radial velocity of 
v t~ o' v {v t -\-cr v ) in the first (fourth) quadrant. The velocity 
dispersion has been estimated to be about 9 km/s for the 
first quadrant and 9.2 km/s for the fourth quadrant [21], 
however it is larger close to the Galactic center. 

The hydrogen number density at tangent points along 
different directions is shown versus height in fig.Q. The 
mid-plane, where the gas density is maximum, almost 
coincides with the Galactic plane. The vertical distribu- 
tion of the gas around the mid-plane can be estimated 
by a Gaussian function of the form exp[— (^^-) 2 ] where 
zo and a z are, respectively, the mid-plane displacement 
and the scale height. The scale height at tangent points 
varies in a small range between 0.1 to 0.2 kpc. 
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FIG. 1: Atomic hydrogen number density at tangent points 
along different directions is plotted versus z. The related 
Galactocentric radius is R = Rq \ sin Z| 



B. Inner Galaxy 

Inside the solar orbit, for a given radial velocity there 
are two heliocentric distances. In order to distribute the 
signal among these points, we assume that the vertical 
extension of the gas at every given R is a Gaussian func- 
tion whose scale height and mid-plane displacement are 
obtained from tangent point with the same R. The point 
whose height is closer to the mid-plane receives more con- 
tribution from the signal and has greater number density. 
The number density at r^, where i indicates either the 
near or the far point, is obtained by 
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C. Local Gas 

At every direction the radial velocities associated with 
negative angular velocities or distances far from the 
Galactic plane (\z\ ^> a z ) are due to peculiar motions 
of the local gas. We assume that the gas with pecu- 
liar radial velocity is locally distributed by a Gaussian 
function with radial scale a r of less than about 4 kpc as 
follows 
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FIG. 2: Atomic hydrogen density map on the Galactic plane. 



where the integral is performed over peculiar velocities 
along the related line of sight. The impact of this as- 
sumption is marginal, since the amount of gas with pecu- 
liar velocity is only 0.033% of the total amount of atomic 
hydrogen gas in the Galaxy. 



IV. RESULTS 

We follow the method described in section Hill to derive 
the number density of atomic hydrogen gas as a function 
of r, I and b centered at the Sun. It is transformed into 
Cartesian coordinates centered at the Galactic center in 
which the Sun is assumed to be at (x, y, z) = (—8.5, 0, 0) 
kpc. The map of number density on the Galactic plane 
z = is shown in fig.(|2|. The distribution of gas on the 
Galactic plane is north-south asymmetric. The density 
peaks at the Galactic center however there is a distinct 
hole right below it. Most of the gas on the Galactic 
plane is concentrated within the radius of about 10 kpc 
whereafter it rapidly dilutes away. 

The global properties of the atomic hydrogen gas dis- 
tribution in the Milky Way are explained in the following 
sections. 



A. The Warp 

There is a large scale warp in the gas disk of the Milky 
Way (see [19] and references therein). The map of the 
mid-plane displacement is shown in fig. Inside so- 
lar circle the mid-plane and the Galactic plane coincide 
pretty well however outside this region the mid-plane is 
warping. The warp is weak up to radii of about 13 kpc 
and then it quickly becomes strong. The mid-plane bends 
up to a height greater than zq = 2.5 kpc at R=20 kpc in 



FIG. 3: The warped Galactic plane. The mid plane displace- 
ment has greater positive (negative) values in darker (lighter) 
regions. 



the north and bends down to zq = —1.5 kpc in the south. 

To better illustrate the warping feature of the gas dis- 
tribution, in figj4]we display the gas number density maps 
at different heights. At negative values of z the dense re- 
gions are in the south and at positive values of z they are 
in the north. 

In fig. ([5]) the mid-plane displacement for different val- 
ues of x is shown against y. The general behavior of in- 
creasing the mid-plane distance from the Galactic plane 
outward the Galaxy is asymmetric with more vertical ex- 
tension in the north. 



B. The Flare 

The balance of gravitational force against the pressure 
force determines the thickness of the atomic hydrogen 
disk. The average scale height, which is defined as the 
distance over which the number density decreases by a 
factor of e, shows a clear flaring [T7[ [22H24] . It increases 
from about 0.2 kpc in the inner Galaxy up to 0.75 kpc at 



R=20 kpc as shown in fig. (16 



top panel). We also show 



the mid-plane density in fig.(6j bottom panel). It peaks 
at the Galactic center and has fluctuations in the inner 
part, then falls down in the outer Galaxy. 



C. Spiral Structure 

The spiral structure can be traced in the surface den- 
sity distribution as regions with over densities [25], [26] . 
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FIG. 4: Atomic hydrogen density maps at different values of 
z. Plots from left to right are, respectively, due to z=-l,l kpc 
on the top panel and due to z=-0.5,0.5 kpc on the bottom 
panel. 
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FIG. 5: The height from the Galactic plane where the number 
density of the atomic hydrogen gas is maximum. 

In fig. ([7]) we show the surface density map, 

E(x,y) = J dzn H (x,y,z) (21) 

in which several spiral arms are evident. There is one 
large spiral arm in the north, the so-called Outer arm. 
In the southern half, the so-called Sagittarius-Carina 
arm close to the solar circle is prominent. The so-called 



FIG. 6: Top: The average scale height versus Galactocentric 
radius. Bottom: The mid-plane density versus R. 

Perseus arm in the south extends to the north and con- 
nects to the Outer arm. 



V. CONCLUSIONS 

In this paper we present a new model for three dimen- 
sional distribution of atomic hydrogen gas in the Milky 
Way. The most recent data on 21cm line emission pro- 
vided by the LAB survey is used. To convert the observed 
brightness temperature distribution to the volume den- 
sity distribution, we assume a purely circular rotation 
curve. We expect it to be a reasonable approximation, 
except in the central region of the Galaxy where the cen- 
tral bar exists. The vertical distribution of gas around 
the Galactic plane is estimated as a Gaussian function. 

The overall structure of the gas distribution discloses 
the warping of the Galactic plane outside the solar cir- 
cle. The bending becomes strong at Galactocentric radii 
greater than about 13 kpc. The thickness of the gas disk 
flares outward the Galaxy. At the same time, the mid- 
plane density falls down in the outer part of the Galaxy. 
Several spiral arms can be traced in the surface density 
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map. It is found that the total mass within a radius of 
20 kpc is 4.3 x 1O 9 M and only 0.033% of that is due to 
local gas with peculiar velocities. 

The derived distribution of atomic hydrogen gas can be 
used to study the propagation of cosmic rays within the 
Galaxy. It can also be used to evaluate the diffuse gamma 
ray maps. Indeed, the structures of the gas distribution 
can be identified by high angular resolution maps of the 
Fermi gamma ray telescope. 

This model is publicly available at the following link 
http : //people . sissa. it/^tavakoli 
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